Spontaneous curvature-induced pearling instability 
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We investigate the instability of a tubular fluid membranes made of a water soluble surfactant. The tubules 
are obtained at high brine salinities. The instability is due to introduction within the vesicle multilayer 
an alkane (Dodecane) . We measure the wavelength of this instability versus the unperturbed radius of the 
tubules. To interpret this dependance we use a model that include not only the surface tension in the elastic 
energy, but the spontaneous curvature as well. The spontaneous curvature is induced by the presence of 
the oil in the bilayer of the membranes, and give a selection of a non-zero wavelength. 

PACS numbers: 68.10.-m, 02.40.-k, 47.20.-k, 87.22.Bt, 
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I. INTRODUCTION 

Vesicles naturally exist in different shapes. The most 
common shapes are spherical and oblates like vesicles, 
toroidal vesicles also exist under some particular condi- 
tions (l]J^l. Tubular vesicles are also common in nature, 
but they are not an equilibrium solution for the elastic 
energy S. One of the most fascinating phenomena are 
the shape changes of cells and the manifold of local insta- 
bilities of the cell plasma membrane in particular and the 
vesicles membranes in general. Under physiological con- 
dition the cell (vesicle) exhibits the familiar biconcave 
shape (the discoide form). Removal of some particular 
proteins leads to the echinocyte form, or to the stomato- 
cyte form The interest of vesicles as model system 
of living cells, for drug delivery or for different industrial 
usage... etc. is the origin of several fundamental studies 
over a couple of decades Local instabilities play an 
essential role for material transport process between ce- 
llular compartments or through the plasma membrane. 
A most important example is the transport of newly syn- 
thesized membranes to the cell envelope, and it involves 
three steps: a) budding of vesicles, b) their fission from 
a parent membrane, c) their fusion with the target com- 
partment. Tubular vesicles when excited by mean of dif- 
ferent techniques (UV radiation, optical tweezers...), ex- 
hibit curious behavior similar to the Rayleigh instability 
in liquid column. This behavior is known as "pearling" 
|@-p[. After an exposition to radiation, tubules become a 
succession of spheres, separated by narrow necks. In the 
case of UV radiation, a long time radiation induces a sep- 
aration of the spheres 0. When the tubules are excited 
by optical tweezers, the shape of the spheres depends on 
the intensity of the laser and a long time exposition of 
the optical tweezers on the tubules induces a separation 
of the spheres separated by narrow necks || . Also it was 
found that vesicles can suffer the same pearling instabil- 
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ity when subjected to a change of the water Ph between 
the inside and the outside of the vesicle || . 

In this article we study an instability that suffer tubu- 
lar vesicles made of a water soluble amphiphile. This 
instability is similar to the ones cited above. However, 
the tubules are not excited via an irradiation, but the 
interlayer distance may be affected by incorporating a 
solvent of the aliphatic chaines of the amphiphiles into 
the layer. In Section II we will describe the experiment 
and in section III we present the theoretical model of the 
instability observed and finally a conclusion is presented. 



II. EXPERIMENT 

It is known that Alkanes are more or less good solvents, 
depending on their chain length, for the aliphatic chain 
of the surfactant called AOT (sodium di-2-ethylhexyl- 
sulfosuccinate) at the oil- water interface AOT is a 
water-soluble surfactant, it forms micelles when dissolved 
in water up to a concentration equal to 30 mmol/L. Be- 
yond this "solubility limit" , AOT forms spherical and 
multilamellar vesicles [jl2| . These vesicles are found to be 
opaque when observed through a microscope and present 
different layer when observed under a phase contrast mi- 
cropscope. We believe they are onions like vesciles. How- 
ever at low salinities and at low AOT concentrations, it 
forms unilamelar sphericals and prolate vesicles as well 
as dumbell shape vesicles. At high salinities (we investi- 
gated the range 0.1-0.175 mol/L), we found that AOT in 
brine forms tubular vesicles |12|. These tubular vesicles 
can be very long and can be of several order of magni- 
tudes of their diameter. There are two ways of preparing 
this phase: Either solubilising AOT in a very short alkane 
(Heptane) , and drying up the solution before adding wa- 
ter on the residue left after the evaporation of Heptane; 
or mixing first AOT and salt in a ration 1:3 in weight, 
then pouring gently the appropriate quantity of water to 
the AOT/Salt mixture to have the desired (AOT/Salt) 
concentrations. For the purpose of this work we used the 
second procedure, so we can avoid any kind of mixing of 
the two alkanes. If we prepare a low salinity solution by 
dissolving the amount of AOT in water and adding brine 
to bring the solution to the desired salinity, no vesicle is 
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observed under the phase contrast microscope. The for- 
mation of these vesicles is probably due to an electrostatic 
effect between the amphiphiles heads [jisfl . In this article 
we will not describe the phase observed when changing 
salinities or amphiphiles concentration. Instead we will 
focus on a particular effect of dodecane on these tubules 
obtained at high salinities. An instability is studied and 
the wavelength selected at the onset of the instability is 
measured as a function of the radius of the unperturbed 
tubules. The selection is found to be a result of the in- 
corporation of the oil into the layer of the tubule which 
causes the spontaneous curvature to deviate from zero. 

The vesicles are observed through an inverted phase 
contrast microscope (40 x objective, Nikon diaphot 200). 
The cell where the vesicles are observed is a 1-mm- thick 
glass cell, the gap between these two glass slides is also 1- 
mm- thick. The AOT-brine solution is left for a couple of 
days, so that any shape due to macroscopique flow -after 
transfering the solution from a test tube to the obser- 
vation cell- desapears, though no shape transformation 
was observed after the transfer of the vesicle dispersion 
into the observation cell. The brine salinities used corre- 
sponding to the minimun of oil-water interfacial tensions 
and are 0.175 mol/L for dodecane and 0.075 mol/L for 
decane and an AOT concentration of 7.5 mmol/L and 4 
mmol/L respectively. The CMC of the same system with- 
out salt is around 2.5 mmol/L. At 0.075 mol/L of NaCl, 
the stable shapes are essentially prolate and spheres; and 
at 0.175 mol/L of NaCl, the shapes are cylinders. This 
shape transformation from low salinities to high salinities 
is probably du to an electrostatic effect, leading to shape 
transformation) |l3 14 1 . We observed the effect of dode- 
cane on stable tubules at 0.175 mol/L (at this salinity the 
Dodecane-Brine interface tension is minimum) ||To[| . The 
introduction of the oil into the cell does not perturbe hy- 
drodynamically the solution and it is a noninvasive way 
to introduce the drop into the solution. To avoids in- 
stabilities such as Marangoni effect, the oil is introduced 
through a less than 1/im diameter orifice. The drop, 
with a volume of the order of 0.25/iL, travells through 
the glasse by capillarity before it reaches the solution. 
The orifice is situated at the lateral wall of the cell, that 
is between the lower glass slide and the piece separating 
the two slides. The drop is then directly introduced into 
the solution. Due to the weak solubility of dodecane in 
water, and its slow diffusion, it takes several hours for 
the instability to start. Most probably the oil travells in- 
side the solution after being incorporated in the swollen 
micelles. After the oil reaches the tubules situated far 
enough from the point where the drop was introduced, 
tubules become unstable, by forming pearls similar to 
the ones observed in Ref . |||| . We have not noticed any 
shape changes when the gap between the horizontal walls 
of the observation cell is much less than 1 mm. Up to 
now, we are not able to know why this effect depends 
on the gap between the cells- walls where the vesicles are 
being observed. A sinusoidal instability is initiated, and 
develops to a peristaltic state, with a reduction in fluctu- 



ations. However large cylinders are stables, this will be 
clarified in the following model. The structure observed 
is periodic and typical necks between the pearls are ap- 
parent. In the case of thick walled vesicles, the pearls 
don't disconnect. However tubules with thin layers, cut 
into separated spheres at the end of the instability. The 
time over which the spheres diconnect is still unknown. 
This kind of state is shown in (fig. [I]). 




FIG. 1. A peristaltic state, in a thick cylinder. The bar 
represents 10 fim. 

Bar-Ziv and Moses |J suggested, after an experiment 
where a tubular phospholipidic membrane was excited 
using an optical tweezers, that, the instability appears 
because of a change in the surface tension induced by 
the tweezers, by analogy to the Rayleigh instability of a 
column of liquid p5|-|l7| . The surface energy should be 
the source of instability, however as the surface energy in- 
creases the linear analysis show that the instability devel- 
ops only longwave modulations without any wavelength 
selection. We will show in the following that in the case 
of our experiment surface tension is not sufficient to se- 
lect a well-defined periodic spatial modulation of tubular 
membranes. Instead, in our case an extra length scale 
becomes important in the problem. We introduce a sim- 
ple model which introduce the physical meaning of this 
length scale. 

Figure (||) , shows the dependance of the dimensionless 
number (q c Ro) as a function of the initial radius of the 
tubule Rq- The wavelength length at the onset hence q c , 
is measured whenever a tubule starts to suffers the insta- 
bility. The dispersion in the data in figure |^ is probably 
due to the fact that the wavelength is not measured ex- 
actly at the onset with an error of the order of 2 seconds 
(that is after the establishement of the perstaltic state). 
The oil molecules reaches the tubules in a random way, 
so the onset is also random. The finer tubules start to 
suffer the instability earlier than the larger ones, this fact 
explains also the dispersion in the data of figure |^. 

The characteristic slope of the experimental values of 
figure |^, represents an extra length (of the order of 2/zm) 
in the problem. The inverse of this length may be in- 
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FIG. 2. The product (q c Ro) as a function of the unper- 
turbed radius (Ro, in microns). The line represents a fit with 
the theory, see later. 



terpreted as a spontaneous curvature c , measuring a 
symmetry breaking between the two sides of the mem- 
brane by incorporating different amount of oil molecules 
(in some sense is the analog of chirality in liquid crystals) . 



III. THEORY 

One may models this by regarding a membrane layer 
as composed of two different layers separated by a dis- 
tance h and with different values of surface tension and 
surface elements one with a + and dS+the other <j_ 
and dS-. Therefore one has an energy of the form 
cr_|_ J dS + + er_ J dS- . This energy would include all 
the elastic contributions (curvature and surface tension 
of the whole layer) |jl8|| . Expanding now both surface 
elements dS± in powers of h, with dS± 
dPdPil + h/2R 1 )(l + h/2R 2 ), one gets 



dl\dl 2 ± 



a.) JdS + (<r+ - a_)h J(±- + ±-\ dS. 



Here ^ 



dVdV 



-k- is the mean curvature and dS 
a geometrical neutral surface element between the two 
layers. Finally, adding a "bending" energy of a tube, the 
free energy is H : 



dS. 



(i) 



Where n is the rigidity, a = (a + + er_) — kcq/2 the "ef- 
fective" surface tension of the bilayer and cq = (cr+ — 



tj-)h/4:K is the spontaneous curvature of the layer and 
is null when the bilayer is symmetric, or when the neu- 
tral surface of the layer is in the middle of the layer [jl9| . 
Here the spontaneous curvature is not introduced as a 
parameter, but as an assymetry in the layer and a shift 
of the neutral surface from the middle of the layer. Also 
the surface tension is expressed as combination of the 
different mechanical parameters of the mambrane.For an 
axisymmetric surface characterized by a local radius R(z) 
the mean curvature is 
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on the 



other hand the surface element 



dS 



R' 2 dz. Here R' = dR/dz 
Let Ro be the radius of a non perturbed tube; as in 
15| we put 



R{z) = Ro 



(3) 



(the square root appears because of conservation of the 
volume of the tube, note that u q= o cannot hold a value 
diffrcnt from zero) into the free energy getting a power 
expansion in u q : 



F = Fo 



k(<?)KI 



(4) 



here 
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Rn 



((1 - c i?o) 2 + 25), 



is a dimensionless number, and the coefficient of 



the quadratic term is 

«(«) = m^(-( c oRo) 2 - (qRo) 2 - <±(c Ro)(qRo) 2 

+ (coRo) 2 (qRo) 2 + 2(qR Q ) 4 
-2S(coRo) 2 (l-(qRo) 2 ))- 



(5) 



The linear instability appears for a(q) < 0. Geometri- 
cally, the transition a(q) = allows us to express the 
control parameter S as a function of the pertubation 
wavenumber qRo- For < coRo < 1/2, S has a minima 
at q = 0. The cylinder is stable for low values of S, how- 
ever, as one increases S up to a critical value (S(q = 0)), 
the tube suffers a long wavelength instability, without 
wave number selection. For 1/2 < coRo < I, S has a 
minima for q = q c . A short wavelength instability ap- 
pears as soon as S > S c see figure (||). From the expre- 
sion (|[) one get explicitly the threshold of the instability, 
hence the wave vector of the most instable mode (q c ), 
and the value of the controle parameter (S) at the onset 
of the instability. 

In the figure (0) , one can see that the control parameter 
vanishes as soon as the radius of the tube is equal to the 
spontaneous curvature, i.e. when the tubules layer are 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 tension vanishes as it should be. 



FIG. 3. The order parameter S as a function of the product 
qRo, the figure in the left is obtained for a zero spontaneous 
curvature, and the one in the right is for Co 7^ 0. 

symmetrical. From ^| and imposing a(q) — dS/dq = 
one gets the critical wavenumber and the threshold. 
Being the first 



q c R = V 1 ± vVl - c R . 



(6) 



Note that there are two different solutions, however the 
sign + represents a situation where q c Ro > 1 and this in- 
stability is possible only if S takes negative valuesQ how- 
ever, the selected wavenumber in this case has a negative 
"mean slope", that is, is a line essentially perpendicular 
to the experimental line given in figure ^. The critical 
surface tension for the case q c Ro < 1 is 

c _ -3 + 2 5 /Vl - c i?o + 4c i? - (c i?o) 2 

The fit represented in figure ^| was done using the ex- 
pression (||), i.e. we tried q c Ro ~ yl — K\/l — cqRq, 
the fit gives K « 1.2 (instead of 1.4) and Cq = 0.16/zra _1 
this value is in the range of value measured for vesicles. 
In figure (||) the value of S c at the q c is of the order of 2. 

At the critical point q c , the sign of the coefficient 
f3(q c ,S c ) of the fourth power term \u q \ , is positive for 
all values of C0R0 such that 1/2 < coi?o < !• There- 
fore this peristaltic instability is supercritical (a second 
order transition). It is interesting to note that for large 
tubular radius there is no instability, as we can see in 
the expression ( |^) where q c becomes a complex number. 
This characteristic is indeed observed in the experiment 
|tt2|| . On the other hand, for small tubules radius, the 



1 Note that this negative surface energy makes instability 
even for Co = 0. 



instability grows for q c = 0, however this is not possible 
because of the conservation of volume, in other words a 
homogeneous growth of the u q is forbidden. In conclu- 
sion in both cases: i) for small tubular radius: Ro < 1/co 
(experimentally we get i?o < 1/im) and for ii) large tubu- 
lar radius Rq > l/c (experimentally R Q > 10/^m) the 
tubules are stables. 



IV. CONCLUSION 

The oil when incorporated in the membrane of the 
vesicle, is distributed assymtrically between the differ- 
ent layers composing the membrane, i.e. between the 
outer layer and the inner layer, giving rise to a non-zero 
spontaneous curvature. Also a possible scenario would 
be that when the swollen micelles reach the tubular vesi- 
cles, they load their amphiphiles at the outer surface and 
a symetry bearking occures between the outer and the in- 
ner layers. A microscopic description of the distribution 
of molecules in terms of the volume fraction of oil intro- 
duced into the alipahtic part of the layer, is necessary to 
understand this effect on the spontaneous curvature of 
the tubules One can wonder of how it is possible 

to vary the spontaneous curvature, without changing the 
length of the surfactant molecule? As for AOT mono- 
layers at the planar interface between the oil and brine 
pha ses, this can be done by changing the brine salinity 
fl22fl , but in the range of the salinity one can probe, the 
Debye length changes from 0.456 nm to 0.608 nm, which 
is not sufficient to change the spontaneous curvature of 
the tubular vesicles radius in the range 5 to 20 /im. This 
can be perfomed using optical tweezers to bring proteins 
to a particular region of the membrane. 

In conclusion, tubules can be excited by different man- 
ners: from UV radiations to laser tweezers; the latter rep- 
resents a powerfull tool to control the onset of the insta- 
bility and the size of the pearls nucleated ||. However, 
it is found that the incorporation of another molecule 
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in this case an alkane in the bilayer, a more or less 
good solvent for the aliphatic chain of the surfactant, 
the tubules suffer an instability similar to the raleigh in- 
stability. Therefore, the non-zero spontaneous curvature 
is created by the incorporation of oil molecules (dode- 
cane) into the membrane of vesicles, inducing a short 
wavelength instability to appear. In our model a surface 
tension of the whole layer is necessary but not sufficient 
to describe the peristaltic states in fluid tubular mem- 
branes. However it is important to introduce an effective 
surface tension (which is probably created by the incor- 
poration of oil molecules between the aliphatic tails of the 
surfactant) composed of the true surface tension of the 
outer and inner layer which suffer elongation and com- 
pression respectively and the rigidity of the layer as well 
as its thickness. 
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